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Abstract 

We compute the fermionic («y) contributions to the flavour non-singlet structure functions in un- 
polarized electromagnetic deep-inelastic scattering at third order of massless perturbative QCD. 
Complete results are presented for the corresponding ny-parts of the three-loop anomalous dimen- 
sion and the three-loop coefficient functions for the structure functions F2 and F^. Our results 
agree with all partial and approximate results available in the literature. The present calculation 
also facilitates a complete determination of the threshold-resummation parameters B2 and D® ls of 
which only the sum was known so far, thus completing the information required for the next-to- 
next-to-leading logarithmic resummation. We find that D® ls vanishes in the MS scheme. 
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1 Introduction 



Structure functions in deep-inelastic scattering (DIS) form the backbone of our knowledge of the 
proton's parton densities, which are indispensable for analyses of hard scattering processes at 
proton-(anti-)proton colliders like the Tevatron and the future LHC. Structure functions are 
also among the quantities best suited for precisely measuring the strong coupling constant a s . 
Over the past twenty years DIS experiments have proceeded to a high (few-percent) accuracy 
and a wide kinematic coverage [Jl|]. More results, especially at high scales Q 2 , can be expected 
from the forthcoming high-luminosity phase of the electron-proton collider HERA at DESY. On 
the theoretical side, at least the next-to-next-to-leading order (NNLO) corrections of perturbative 
QCD need to be taken into account in order to make full use of these measurements and to make 
quantitatively reliable predictions for hard processes at hadron colliders. 

The calculation of NNLO processes in perturbative QCD is far from easy. For deep-inelastic 
structure functions, in particular, the current situation is that, while the coefficient functions are 
known to two loops [EL BL HL B| H], only six/seven integer Mellin moments of the corresponding 
three-loop anomalous dimensions have been computed for lepton-hadron [|7|- D- 0] and lepton- 
photon DIS [pH, together with the same moments of the three-loop coefficient functions. The 
hadronic results have been employed, directly [ |TT| , [12| , [13], [14]] and indirectly [[15|, [16]] via x- space 



approximations constructed from them [jT7|, [T8], [T^], to improve the data analysis and some hadron- 
collider predictions. However, the number of available moments is rather limited, and hence these 
results cannot provide sufficient information at small values of the Bjorken variable x. 

For the complete information one needs to obtain either all even or odd (depending on the quan- 
tity under consideration) Mellin moments, or do the complete calculation in Bjorken-x space. We 
have adopted the first approach. Following the formalism of ref. [ |20| , [21] , 0, [|] to obtain the lower 



fixed moments, we have used recursive methods to extend the calculation to all values of the Mellin 
moment N. This is by no means trivial, since before the start of the calculation the mathematics 



of the answer was still poorly understood Q22|]. Hence it was first necessary to develop the under- 
standing of harmonic sums [23], |25|, |2^] and harmonic polylogarithms [ LT7| , [2~8| , ^9p. In addition 
the Mellin transforms and the inverse Mellin transforms from Bjorken-jc-space to Mellin space and 
back had to be solved [ ]2"9| ] . These conceptual problems have been overcome and the method has 



been shown to work in a complete re-calculation of the two-loop coefficient functions [|5p . 

The concept of working in Mellin space is not new. This method was already used in the early 



QCD papers Q30| , |31] , 32 ], But even in the case of the two-loop anomalous dimensions it was still 



possible to do the resulting sums in a rather direct manner ^3p. This changed with the two- 



loop evaluation of <5i/<5t in which Kazakov and Kotikov 034] , p5| ] managed to obtain some of the 



integrals via recursion relations or first order difference equations. 

In this article, we present the fermionic (ny) corrections to the flavour non-singlet structure 
functions in electromagnetic DIS at the three-loop level. This includes the three-loop anomalous 
dimensions which are needed for the completion of the NNLO calculation, as well as the three- 
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loop coefficient functions which are, at least at large x, the most important contribution to the 
next-to-next-to-next-to-leading order (N 3 LO) correction [|36|]. Of course, the rc^-part is not the 
complete calculation. Yet we decided to present it already now, because the complete calculation 
(including the singlet part) will still take quite some (computer) time. Thus, for the time being, 
NNLO calculations of the Drell-Yan process [TJ7|, [38|] (which are relevant for luminosity monitoring 
at Tevatron and LHC ||3^, |Mj[ [41]]) and of Higgs production [|38|, |2|] have to rely on parton 
distributions evolved with the approximate splitting functions of ref. [JT^]. Our present calculation 
provides a first check of the reliability of these approximations. More importantly, it turns out 
that already this calculation is sufficient to provide the last relevant missing information for the 
extension of the soft-gluon (threshold) resummation for DIS [B3L B3L 03p to the next-to-next-to- 



leading logarithmic accuracy ffifa. In fact, our result for the resummation parameter D® ls is most 
intriguing and calls for further studies. 

This article is organized as follows. In section ^ we outline those parts of the calculation, which 
differ from previous two-loop [^p and fixed-Af three-loop [[7], [5p analyses. The present calculation 
does not yet involve the full complexity of the method, as the most difficult diagram topologies do 
not occur. Therefore we postpone a full account to a later publication. In section |3| we present our 
explicit even-Af Mellin- space results, except for the rather lengthy expressions for the three-loop 
coefficient functions which are deferred to appendix |A] The corresponding three-loop quantities 
in jc-space can be found in section [|. Here, instead of writing down the cumbersome exact ex- 
pressions for the coefficient functions, we follow the procedure applied in refs. [ [T7| , [T8|] to the 
two-loop coefficient functions, and provide approximate parametrizations which are compact and 
sufficiently accurate for all numerical applications. In section |]we then discuss the implications 
of our calculation on the soft-gluon resummation, before we summarize our results in section ^|. 



2 Method 

Because we are considering only the non-singlet structure functions in this paper, the method for 
the calculation of their moments closely follows ref. [0]. Hence there is not much need to explain 
the physics of the method here again. Thus we will discuss only the differences introduced by the 
fact that we now compute all moments simultaneously as a function of the moment number N. 
Since N is not a fixed constant, we cannot resort to the techniques of ref. [|7|], where the Mincer 
program [^7|, ^Bp was used as the tool to solve the integrals. Instead, we will have to introduce new 



techniques. However, we can give N a positive integer value at any point of the derivations and 
calculations, after which the Mincer program can be invoked to verify that the results are correct. 
From a practical point of view this is the most powerful feature of the Mellin-space approach, as it 
greatly simplifies the checking of all programs. 

Similar to the fixed-A^ computations of refs. 0], the diagrams are generated automatically 
with a special version of the diagram generator QGRAF p9p . For all the symbolic manipulations 



of the formulae we use the latest version of the program FORM [pOQ. The calculation is performed 



2 



in dimensional regularization [ pl| , |52| , p3[ p4| ] with D = 4 — 2e. Hence the unrenormalized Mellin- 
space results will be functions of £, N, and the values £3 5 of the Riemann ^-function. The 
renormalization is carried out in the MS-scheme [|55], |5§] as described in ref. |[7p. 

We distinguish three categories of diagrams: complete diagrams, composite building blocks 
and basic building blocks. A complete diagram is a Feynman diagram with all its structure like 
traces and dotproducts in the numerator. Such a diagram may lead to a large number of more 
fundamental integrals that cannot be reduced by considerations of momentum conservation only. 
For the understanding of composite and basic building blocks, one has to realize that in the frame- 
work of the operator-product expansion we eventually have to take N derivatives with respect to 
the parton momentum P after which P is put equal to zero. This projects out the iV-th Mellin 
moment [57p and it effectively amputates the legs of the parton, leaving us with propagator-type 



diagrams. Therefore, we define the topology of a diagram as the propagator topology when the 
P-momentum legs have been amputated. The three-loop propagator topologies of the BE (Benz) 
type and of the 04 type are shown in fig. [T]. For the notation we refer to refs. Q4"7|, |4"%|] . The external 
lines in the propagator topology are referred to as Q. 




Figure 1: The topologies BE (left) and 04 (right) of propagator-type diagrams, with the line num- 
bering as employed in figs. ^ and [3| below. The external lines carry the momentum Q. 



When the numbers of the position of the incoming and outgoing momenta have been attached 
we are referring to subtopologies. For instance, BE13 is a subtopology of type BE in which the 
momentum P comes in in line 1 and leaves in line 3, assuming the numbering of the BE topology 
as in fig. []]. We define basic building blocks (BBB) as integrals in which both the incoming and the 
outgoing P-momentum are attached to the same line as, e.g., in BE22- In composite building blocks 
(CBB), on the other hand, the incoming and the outgoing P-momentum are attached to different 
lines, as in the case of BE13 mentioned above. Of course we could have introduced names for all 
these three-loop four-point functions, but since eventually the P-momentum legs get amputated the 
above notation seems the clearest scheme. In this way, it is only a small step to an easy pictorial 
representation of the integrals as used in ref. [|p. 

For the calculation of the n/-parts of the non-singlet structure functions F2 and Fl, one does 
not need to consider all three-loop topologies. In fact, the only genuine three-loop subtopologies 
one has to solve are of the BE-type and of the 04-type, with two complete diagrams of either type 
entering the calculation. These diagrams are displayed in figs. [2| and 
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Figure 2: The diagrams of topology BE which contribute to the fermionic part of the non-singlet 
structure functions F2 and F^. The subtopologies are BEiq (left) and BE13 (right). 




Figure 3: The diagrams of topology 04, which contribute to the fermionic part of the non-singlet 
structure functions F2 and Fl. These diagrams are examples of the subtopologies 04iq and 04ig. 



The main problem we are faced with as compared to the corresponding two-loop calculation [§] 
is that the necessary reduction equations become much more complicated. In addition the set of 
equations that was available at that moment was not maximal, and the structure of the more com- 
plicated topologies needs the maximal set of reduction equations. This involves some equations 
in which the vanishing of P ■ P is an issue that should be considered with care, a point that will be 
explained in full detail in a later publication in which we have to deal with all topologies. 

For this calculation we first studied the basic building blocks. Here it is easy to expand the sin- 
gle propagator that contains the momentum P to sufficient powers in P for the N-th moment. This 
number of powers can be less than N as there may be some powers of P in the numerator already. 
It should be noted that if we have a power of P ■ Q in the numerator we are effectively computing 
the moment N— 1 of the integral. At this point we write down all equations based on integration 
by parts, all scaling relations and all form-factor equations that can be constructed. Next follows 
a potentially difficult process in which we have to combine these relations to construct equations 
that can systematically bring the powers of the denominators in the integral down, either reducing 
them to zero or leaving them at a fixed unique value. When a line is eliminated a simpler topology 
is reached and we can refer to the reduction equations for that topology. This will eventually lead 
to a topology that is simple enough so that we can calculate the integral. 

A problem arises when a reduction equation cannot take the power n of a denominator beyond 
one, because it is of the type nl(n + 1) — >• /(ft). Also when the power of a propagator is not 
an integer, one can only try to reduce this power to a fixed value for which the line does not 
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vanish. In these cases we leave the power either at one or at this fixed value (usually 1+8) and we 
continue with the next propagator. Eventually we will have integrals in which only the line with 
the power that involves N will not have a standard value. The N will occur both in the power of 
the denominator p ■ p and in the numerator IP ■ p when p is the momentum associated to this line. 
By this time our remaining equations which have also been treated to bring the powers of their 
progagators to standard values may have become rather lengthy. But we still need them to bring 
the difference of the powers in the numerator and the denominator to a fixed value, which can be 
either 0, 1 or 2, something we have to leave open for reasons explained below. Now there are 
several possibilities: If there is only one integral left of the type to be solved, the equation directly 
determines the solution. This is however rarely the case. Usually there are several terms left, each 
with a different power of IP ■ Q, leaving us with an equation of the type 

ao(N)I(N)+ ai (N)I(N-l) + ... + a m (N)I(N-m) = G(N) (1) 

in which the function G refers to a potentially horrendous combination of integrals of simpler 
topologies. Eq. (|T]) defines a recursion relation or difference equation of order m. For the present 
calculation we did not have to go beyond order 2. It should be noted that recently difference 
equations have been encountered by Laporta in refs. []5£ 



These difference equations can be solved by making the ansatz that the solution will be a 
combination of harmonic sums. If the proper combination has been selected, the coefficients of 
the harmonic sums can be obtained by substituting the trial solution into the equation and solving 
the resulting, potentially large set of linear equations. There can be several thousand equations in 
such a system. Of course one has to have a solution for the function G(N) in eq. ([[]) first, and 
m— 1 boundary values are required. Because these boundary values are basically fixed integer-Af 
moments, they can be obtained using the Mincer program. The need for knowing the function G 
puts a rigid hierarchy in the order in which we have to treat the topologies. 

Solving the difference equations is rather slow work. Hence we compute their solutions only 
once and tabulate the results. We usually do this for several values of the difference of the powers 
of the numerator and the denominator as mentioned above eq. ([]]). The reason behind this is that 
the equations we use for either raising or lowering this difference may contain a so-called spurious 
pole in £ when we try to bring this value to one. The concept of spurious poles is rather important. 
The rule of the triangle [150, BTH, for instance, can involve a factor proportional to 1/e when an 



integral is being reduced. Close inspection reveals that when powers of the loop momentum are 
present in the numerator, it is possible that more than one of such poles is generated before a 
denominator is removed. The resolved triangle [ |62| ] shows, however, that it is possible to sum 
all contributions of such a reduction and that eventually there can be no more than one pole per 
eliminated line. This means that the extra poles should cancel between the many generated terms. 
But when we work only a to given cutoff in powers of £ (both for reasons of economy and because 
we cannot evaluate some integrals easily beyond certain powers in £) these temporary poles could 
spoil the final result in the same way as such things can happen in numerical calculations at fixed 
precision. We call them spurious poles, because in principle they can be avoided. One of the 
greatest difficulties in deriving reduction equations is to indeed avoid such spurious poles. For 
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some integrals, no spurious-pole free formula could be found for the last reduction when bringing 
the difference of the power of the denominator and the numerator either from 2 to 1 or from to 1 . 
This problem has been circumvented by solving the resulting difference equations for the three 
cases 0, 1 and 2 separately. 

Since the evaluation of all the basic building block integrals that can occur requires very much 
computer time, we decided to tabulate all these integrals for the complicated topologies. This saves 
much computer time, because each integral is typically used many times. Also the more compli- 
cated topologies use many integrals of a less complicated type, hence their evaluation becomes 
much faster once the latter integrals have been tabulated. 

The next step is the evaluation of the composite building blocks. Here again we first construct 
all possible equations. It turns out to be most economic to leave the propagators with the momen- 
tum P unexpanded. The fact that eventually an expansion to N powers of P will take place then 
requires some special calculational rules. If for instance an equation is multiplied by P ■ Q, we 
have to replace N by N— 1 at the same time. Such rules can be major sources of errors. Hence it is 
very fortunate that at any moment we can decide that N has a fixed value like four or five and then 
evaluate the integrals in the equation with the Mincer program to see whether it is still correct. 

The equations are used to set up a reduction scheme that is similar to the one for the BBB's. If a 
line is eliminated we obtain either a simpler topology or a BBB. Otherwise we reduce the power of 
a denominator to a standard value and continue with the next line. In numerous cases a reduction 
can only be done by means of a difference equation ([]]). If this is a first order difference equation it 
can be solved directly, introducing one sum. Such sums are of a benign type and can be evaluated 
afterwards. If the difference equation is of a higher order we have to consider all more fundamental 
integrals first before we can solve the equation. The further reduction scheme becomes then rather 
complicated, but not impossible. On the average each subtopology can require several weeks of 
work before it has been completely solved by these methods. 

Analogous to the BBB case discussed above, we have tabulated the more complicated CBB 
integrals entering the calculation. It is not only a matter of computer time that renders this neces- 
sary. Also the size of the intermediate expressions becomes a most relevant factor: if one is not 
careful, even a hard disk of 100 GBytes can become restrictive. In practice, already expressions 
significantly larger than 10 GBytes took too long for evaluation without further optimization. A 
careful hierarchy of tabulation managed to avoid these problems. 

Having programs for all basic and composite building blocks renders the remainder of the 
calculation rather straightforward. The major difference to the fixed-moment calculations is now 
that we obtain much longer results due to the presence of the parameter N in the answer. We have 
checked the correctness of each individual diagram for several values of N, by comparing with the 
results of a Mincer calculation. In addition we have compared the complete renormalized results 
with the results in the literature for the available values of N. 
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3 Results in Mellin space 



Here we present the iV-space coefficient functions and the anomalous dimensions up to the third 
order in the renormalized coupling a s . All results are given in the MS-scheme with the renormal- 
ization scales identified with the physical hard scale Q. Thus the perturbative expansion of the 
non-singlet coefficient functions and anomalous dimensions can be written as 

(Wo.,*) = £ (2) 

n—Q 



y m (a s ,N) = £ (%.y + \S(N) (3) 

H=0 



with i = 2,L in eq. (Q). There is no need to consider different choices of the scales, as they do not 
require functions beyond those introduced in eqs. (^|) and (T3|). Recall that for F 2 the (n + l)-loop 
anomalous dimensions and the n-loop coefficient functions together form the N n LO approximation 
of (renormalization-group improved) perturbative QCD. 

There are, actually, three different non- singlet combinations of coefficient functions and split- 
ting functions. These combinations all coincide at order a s , but they all differ beyond the second 
order. Only the so-called '+' -combinations (involving sums over quarks and antiquarks) are probed 
in electromagnetic DIS, hence only these quantities are addressed in the present article. Conse- 
quently our results below apply directly only to all even-integer values of N from which, however, 
the results for arbitrary N can be uniquely inferred by analytic continuation. 

Our Af-space results are expressed in terms of harmonic sums S^(N). In the following all 
harmonic sums are understood to have the argument N, i.e., we employ the short-hand notation 
Sff t = Sff l (N). In addition we use operators N± and N±\ which shift the argument N of a given 
function by ± 1 or a larger integer i, 

N±/(A0 = f(N± 1) , N ±i /(A0 = f(N±i) . (4) 

We normalize the trivial leading-order (LO) coefficient function and recover, of course, the 
well-known result for the LO anomalous dimension [BUI RTJ 



42. (#) = X ' (5) 
Yi?(A0 = C F (2(N_+N+)5!-3) . (6) 



In our notation, the next-to-leading order (NLO) non-singlet coefficient function for F 2 D56Q and 



the corresponding anomalous dimension [ p3j p3| ] read 

4% (#) = c f0 n +^i + 2Si - 9 + (N_ + N+) [25 u - 3Sj - 2S 2 ]), (7) 



y$(N) = AC A C F (2N + 5 3 - — - 25_ 3 - ySj + (N_ +N+; 



S\ -\-1S\ _? Si 

18 ' 6 



7 



+ 4C F n / (^ + ^S 1 -(N_ + N + ) 



T 1 " 3 52 



+ 4C F 2 (45_ 3 + 25i + 25 2 - 



S 2 + 2S 3 



-(N- + N+) 



Si+4S. 2 + 2S,2 + 2S 2 ,+S3 



(8) 



At next-to-next-to-leading order (NNLO, N 2 LO), we have re-calculated the coefficient function 



(2) 

c 2 ns °f re ^ s - IB §P' an ^ we nave com P u ted all fermionic contributions to the splitting function y 
the n\ terms being a new result of this article. The fermionic NNLO corrections — the complete 



(2) 
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expression for the N- space coefficient function can be found in ref. [H] — are given by 



41W =4C F n f ((l-N 



122 7 
L27 5l + 6 5M 



13 



1 



1 



"(N- 
457 



Y^(iV) = 16C A C F / I/ ^3--4-— 5_ 3 -— 5 3 + -5 1) _ 2 --5_4 + 25i i i-— S 2 + 



6 5l ' 1+ 144 
10 



-1) 
247 



89 
108 



Si-S 2 



19 

Si + — N+S 2 
108 6 + 



(N_+N- 



25 



rS 3 



(9) 



257 



3^-3,1 



■N- 



2 2 " 
S 2 1 — —S3 1 — — S4 

' 3 ' 3 



+ (1-N_ 



23 ' 
LIS 53 - 52 



11 317 o 16_ 2 c 1 
+ IS* 3 " 108 52 + ~9 Sl - 2 ~ 3 51 - 2 ' 1 " 3 51 '- 3 



(N_+N 
1 



27 

1237 ^ 
Su + ^r^Si 



Si 



-^1,3 -2^.1 -3 
2 c 11 1 

9 5l "54 52 + T8 53 

4 59 4 20 20 8 8 4 

+ ~ S_ 3 1 - — S 2 + -5_ 4 - ^S_ 3 + —Si - -Si _ 2 - -Si 1 - -Si 2 + N. 
3 ' 36 3 9 9 3 ' 3 ' 3 ' + 



216 

S 2 - 2 + Si^3 + -S3 1 



+ 16C ^/(^-^ 5l + ^ 2 + ^- +N - 



+ (1-N_ 



67 



:S 2 



4 4 



+ 16CA/ (|_| & 

[-25 4 
yS 3 --S 3 ,i 

+ (N_+N + ) S^3-^Si-^Si,_ 3 + fsi,- 2 



.36 

4 r 4 r 16 _ 4 c 11 2 c 10 _ 1_ 2 r 8 r 

- 3*1,-2,1 + 3^1.1 + -g i> 2 ~ 3 1,3 + 18 ~ 3 2 ~ 2 + "9" 2,1 + 2 ~~ 3 2,2 ~ 9 



(10) 



The corresponding formulae for the longitudinal coefficient function Cl.iis are deferred to the 
appendix, together with the rather lengthy /V-space results for the fermionic parts of the third-order 
coefficient functions c- ns , i = 2,L, which partly also represent new results of this article (the nf 

term for F 2 has already been presented in ref. [)6"3"1]). Notice that c L 4 can be considered a NNLO 

quantity, since C F vanishes at order a!?. On the other hand cl 3 ^ represents, at least at large N, the 
dominant part of the N 3 LO corrections to F 2 P^|]. 



As briefly mentioned at the end of section ^, we have subjected our results to a number of 
checks. First of all, we have calculated some lower even moments in an arbitrary covariant gauge 
with the Mincer program [^7], |4H[ ], keeping the gauge parameter £, in the gluon propagator. All 
dependence on £, does cancel in the final results. Secondly the ^^-contribution to y$ is known 
from the work of Gracey [B4| and we agree with his result. Furthermore the coefficients of \n k N, 

(3) ' 

k = 3, . . .5, of c\ 4(^0 agree with the prediction of the soft-gluon resummation [|53|]. Finally, we 
have checked the result of each individual diagram for several integer values of N by comparing 
with the results of a Mincer calculation. Thus, as the strongest check, our results reproduce the 
fixed even moments N = 2, . . . , 14 computed in refs. [|7], [8|, 
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4 Third-order results in x-space 



The x-space coefficient functions and the splitting functions are obtained from the results of the 
previous section by an inverse Mellin transformation, which maps the harmonic sums of moment 



space [ |23| , |24] , [25], £6]] to harmonic polylogarithms in x-space [ [271 , P8| P9Q - This transformation can 
be performed by a completely algebraic procedure Ep based on the fact that the harmonic sums 



also occur as coefficients of the Taylor expansion of harmonic polylogarithms. 



Here we confine ourselves to the third-order results; for the two-loop non- singlet splitting 
functions and coefficient functions the reader is referred to refs. [ |671 0, §]. For brevity the exact 
results are written down only for the splitting function, conventionally related to the anomalous 
dimension (Q) by 



dx x 



N-l nW, 



(id 



>(«) { , 



The fermionic part of i 5 ^ involves only simpler harmonic polylogarithms which can be expressed 
in terms of the usual (poly-)logarithms. The x-space analogue of eq. (|T0|), graphically displayed in 
fig. ffl, can thus be written as 



3„ 167. _ 1. , _ 1. 2 



7 j 6 2 ^ 3 " 108 ln W + 3 Hx) *> 2 ~4 ln W ln( 1 " X) 



1 



1 



1. 



1 



— ln 2 (x) - —In \x) - ~ln(x)Li 2 (x) + -Li 3 (x)j +p qq (-x) [-£ 3 - - -ln(l +x)£ 2 
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+ 7ln(x)^ 2 - ^ln(x)ln(l +x) + — ln 2 (x) - ^ln 2 (x) ln(l +x) + -^-\n 3 (x) - ^Li 2 (-x) 
ri 18 6 i o ri 



1 

18 



10, 

t 



i 

6 

- Vi 3 (-x) - ^Li 3 (x) + ^H_i A i(x)l + (1 +x) + ^ln(x) - ^Li 2 (x) - ^Li 2 (-x) 
J J J J Lb 223 

2 1 i r1 257 17 1 

- -ln(x)ln(l +x) + -ln 2 (x)] + (1 -x) - +ln(l -x) - - ln(x) - -ln 2 (x) 



+ 8(1 -x) 



167, 



1 



25, 



20 



.4 54 

5 . , , 1. , _ 10 



18 



k + ^ + ^Ca +16CW p qq (x) ln(x)-- + -ln 2 (x 



1 



1 



17 



L 144 27 



.54 



-— +l6C F \[ Pqq Lv) -C3 



54 36 

5 55 

3^ 3 48 

2 



12 

+ ^ ln(x) + ^ ln(x) £ 2 + y ln(x) ln( 1 - x) + - In 2 (x) + - In 2 (x) ln( 1 -x) + 1 ln(x)Li 2 (x) 



^Li 3 (x) - ^ln 3 (x)j +p qq (-x) - C 3 + ^ln(l +x)^ 2 - ~ln(x)C 2 - ^hi 2 (x) 



1 

18 



10, 

t 



4 

3 



1 

3 



20 1 1 20 2 2 

+ — ln(x)ln(l +x) - -ln 3 (x) + -ln 2 (x)ln(l +x) + —U 2 (-x) + -Li 3 (-x) + -Li 3 (x) 



7 67 4 12 

— ln 2 (x) - — ln(x) + -ln(x)ln(l +x) + — In 3 (x) + -Li 2 (x) 



-1h_i A i(x)]+(1+x) 
+ 1 Li 2 (-x)] + (1 -x) [I ln(x) - y - 1 ln(l -x) + ~ ln(x) ln(l -x) - I ln 2 (x) 



-8(1 -x) 



23 5 „ 29 r 2 17 
Ll6 " 12^ 2 " 30^ 2 + _ 6^ 3 



(12) 
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where we have introduced 

p qq (x) = 2{\-x)- l -l-x (13) 

and all divergences for x — > 1 are understood in the sense of +-distributions. In eq. dT2"| ) we have 
left one particular harmonic polylogarithm, H_i o,i(jc), unsubstituted. This function is given by 

fX Ay J 

H_i,o,i(*) = / -r— Li 2 (z) = U 2 (x)\n(l+x) + -S L2 (x 2 )-S h 2(-x)-S ia (x) , (14) 
Jo l+z I 

where the representation by the Nielsen function Si ,2 has been derived in ref. [|55p. H-i^ifc) can 
also be expressed in terms of trilogarithms, albeit with more complicated arguments [§]. 



-180 
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(1-x) P + ( j } (X) 




N = 2...12 
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X 



-1 - 
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(l-x)P + f (X) 
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0.2 0.4 0.6 Oi 
X 



Figure 4: The and nj parts \{x) and P\_ 2 (x) of the three-loop non-singlet splitting function 
(|12|), multiplied by (1 — jc) for display purposes. Also shown in the left part (dashed curve) is the 
uncertainty band derived in ref. [JTPp from the lowest six even-integer moments \U\ BJ, H]. 



The x-space coefficient functions involve harmonic polylogarithms of weight four, which in 
general cannot be expressed in terms of standard polylogarithms and Nielsen functions anymore. 
Instead of writing down the cumbersome exact expressions, we prefer to present sufficiently accu- 
rate, compact parametrizations in terms of the +-distributions and end-point logarithms 



ln k (l-x) 



1 



-x 



L\ = ln(l — x) , Lq = \nx 



(15) 



It is convenient to apply this procedure (which has been employed in ref. [ |T7| ] for the two-loop 
coefficient functions) also to the nj part of the splitting function (TT2|). Inserting the numerical 
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value of the QCD colour factors, this function can be approximated by 



P^\x) £ nA- 183. 187 2) - 173.927 8(1 -x) -5120/81 Li - 197.0 + 381. lx + 72.94 x 2 



+ 44.79 x 3 - 1.497x4 - 56.66 LoLi - 152.6 Lo - 2608/81 h\ - 64/27 4) 

7 



+ nh-Vo -(51/16 + 3^3 -5C2) 8(1 -x)+x(1-x)-%(3/2Lq + 5) + 1 



+ (l-x)(6+ll/2L + 3/4L 2 ) ))64/81 . (16) 
Corresponding parametrizations for the three-loop coefficient functions read 

c^ s (x) ^ n f { 640/81 © 4 - 6592/81 ©3 + 220.573 © 2 + 294.906 ©1- 729.359 © 
+ 2572.597 8(1 -x) -640/81/4 + 167. 2L 3 , -315.3L 2 + 4742Li 
+ 762. 1 + 7020x + 989.4x 2 + L L X (326.6 + 65.93 L + 1923 L\ ) 
+ 260. 1 L + 186.5 4 + 12224/243 L 3 + 728/243 Lq ) 
+ n/ (64/81 ©3 - 464/81 © 2 + 7.67505 <D X + 1.00830 £> - 103.2655 8(1 - jc) 

- 64/81 L\ + 15.46 L\ - 51.71 L\ +59.00 x + 70.66x 2 +L Li (-80.05 

- 10.49 Lo + 41.67Li)-8.050Lo-1984/243Lj-368/2434) , (17) 

c 2s W - "/ ( 1024/81 Lf - 1 12.4L 2 + 340.3 L x +409 - 210x - 762.6x 2 - 1792/81x4 
+ L Li (969.2 + 304. 8 L - 288.2 L x ) + 200. 8 L + 64/3 L 2 , + 0.046 8( 1 - x) ) 
+ nj ( 3xL 2 + (6 - 25 x) L x - 19 + (317/6 - 12 £ 2 ) x - 6xL Li + 6xLi 2 (x) 

+ 9xLo-(6-50x)Lo) 64/81 . (18) 

The n 2 parts of Pns and c Lm , the +-distribution contributions (up to a numerical truncation of 
the coefficients involving £,•), and the rational coefficients of the (sub-)leading regular end-point 
terms are exact in eqs. dTBj ) - (JTS|). The remaining coefficients have been determined by fits to the 
exact results, for which we have used the Fortran package of ref. []68[j. The above parametrizations 
deviate from the exact expressions by one part in thousand or less, an accuracy which should be 
amply sufficient for foreseeable numerical applications. 



5 Implications for the threshold resummation 

The large-TV/ large-x behaviour of the three-loop splitting functions and coefficient functions is of 



special interest in connection with the soft-gluon (threshold) exponentiation [|43J [44], [45]] at next-to- 
next-leading logarithmic (NNL) accuracy. Here the coefficient function for F 2 ns can, up to terms 
which vanish for N — > °°, be written as 

C 2 , ns (a s ,/Y) = (l+a s g 01 +a 2 g 02 + ...) exp[Lg 1 (a s L)+g 2 (a s L)+a s g 3 (a s L) + ...] (19) 

with a s = a s /(47t) and L = \nN. The functions g l depend on (universal) coefficients A,</ and 
#(•</_! and process-dependent parameters Df^ {l as described in ref. [|4"6j], where also the explicit 
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expressions for the functions g 1 2 3 can be found. Hence the NNL function g 3 involves the new 
coefficients A3, B2 and D^. These coefficients can be fixed by expanding eq. ( JT9| ) in powers of 
a s and comparing to the result of the full fixed-order calculations. 

In the MS scheme adopted in this article, the parameter A3 is simply the coefficient of In N in 

(2) (2) ' 

Yns (N) or, equivalently, of 1/(1 — x)+ in P^s (x). Its fermionic part is thus known from eq. ([12]), 

112 



A3 



= CaCfu 



836 160 



<3 









16" 


+ Cprif 




+ Cpnj 


~27_ 



(20) 



The numerical value can be read off from eq. (fl6|). Like for the whole of P\J{%), as shown in 
fig. |], this result is consistent with, but supersedes the estimate derived in ref. [ |T9"| ] from the first 
six even-integer moments. Parallel to our work A3 L has also been calculated in ref. [|69|]. 



The combination B2 + D^ 1S has been determined in ref. [|4^] by comparing the expansion of 
eq. ( |I9| ) to the InN term of the two-loop coefficient function of ref. [§]. As the ln 2 ./V (or 
©1) contribution to ci^ s involves a different linear combination, Po(#2 + 2D^ IS ), of the very 
same coefficients, B2 and £>P can be disentangled using the three-loop coefficient function. The 
analytic results for the two new +-distribution coefficients read 



,(3) 
•2,ns 



<D x n f 



Cp^Cpn 



f 



15062 512 „ _ 

1 C2+ I6C3 

81 9 s s 



83 1£0f . 112. 
— + 168^2 + -^^ 



.(3) 
"2,ns 



■D n f 



+ CpHf 

C A C F n f 
4226 



940 



81 

160906 
729 



32 * 
y^ 2 



(21) 



9920 
~8T 



776 „ 208 r 2 



27 



^2-60^ + 16^ 



■f 



15 
8714 



2003 



108 



232 y 
+ — £2 



729 27 



32 
27 



C3 



(22) 



(the coefficients of £> 2 5 can be found in ref. [|65|]). In fact, due to the prefactor f$o, the complete 
results for B2 and can a l rea dy be inferred from fermionic result (|2T|), yielding 



5, 



---24^ + 12^2 



+ CfCa 



3155 ^„ 44/ 
-— + 40^ + ^2 



D 



DIS 



. 



247 8 
27 3^ 2 



(23) 
(24) 



The vanishing of Df ls and DJ? 1S — in contrast to the Drell-Yan process, where D2 is different 
from zero [ ff5H — calls for a deeper explanation, possibly offering an all-order generalization. 

Finally we note that, once the non-fermionic contributions to the 3-loop non-singlet splitting 
functions and coefficient functions are completed, the NNL threshold resummation facilitates a 



prediction of the first six towers of logarithms, i.e., the coefficients of a" In 
C2,ns at all orders n > 3. We will return to this issue in a later publication. 



2n—i 



N, i = 0, .. .,5, of 
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6 Summary 



We have computed the fermionic (n^-enhanced) third-order contributions to the structure functions 
Fi and Fi in electromagnetic deep-inelastic scattering. The calculation has been carried out for all 
even-integer Mellin moments N, by solving the three-loop integrals by means of recursion relations 
(difference equations) in N. This progress with respect to previous computations restricted to some 
fixed moments N is especially due to an improved understanding of the mathematics of harmonic 
sums and difference equations, and the implementation of corresponding tools in the symbolic 
manipulation program FORM which we employed to handle the huge intermediate expressions. 
We are confident that our approach will enable us to compute all three-loop corrections in DIS. 

We have thus been able to derive the complete expressions for the corresponding fly-parts of 
the NNLO anomalous dimensions and splitting functions and the N 3 LO coefficient functions for 
F2 and Fi. The results have been presented in both Mellin-A^ and Bjorken-x space, in the latter case 
we have also provided easy-to-use accurate parametrizations. Our results agree with all partial and 
approximate results available in the literature for these quantities, in particular we reproduce the 
even-integer moments N = 2, . . . , 12 computed before. 

The present results for the three-loop splitting function represent a step towards completing 
the ingredients required for NNLO calculations of hard-scattering processes involving initial-state 
hadrons in perturbative QCD. The three-loop coefficient functions for the most important structure 
function F2 form the dominant part of the N 3 LO corrections at large x, thus facilitating extractions 
of a s with a distinctly reduced theoretical uncertainty. Already the fly-part computed in this article 
leads to a complete determination of the threshold-resummation parameters B2 and D® ls — includ- 
ing the non-fermionic contributions — of which only the sum was known so far, thus practically 
completing the information required for the next-to-next-to-leading logarithmic resummation. 
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A Appendix 



All results for the non-singlet anomalous dimensions and coefficient functions presented in this 
article can be obtained as a FORM file from the preprint server |http: //arXiv.org| by down- 
loading the source file. Furthermore they are available from the authors upon request. 
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The fermionic parts, i.e., all terms proportional to CaCjuij, Cp-n^ and CfHj of the three-loop 
coefficient function for the electromagnetic structure function F2 are given by 
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(A.3) 



The contributions to the 3-loop longitudinal coefficient function corresponding to eq. ( |A.1| ) read 
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Recall that in all our formulae the expansion parameter is normalized as in eq. (@). The operators 
N± and N±, have been defined in eq. @). 
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